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ANALYSIS OF IONOSPHERIC REFRACTION 
ERROR CORRECTIONS FOR GRARR SYSTEMS 

A, J* Malinckrodt, H, C* Parker, J. H* Berbert 


ABSTRACT 

This report presents a determination of the ionospheric refraction correction require- 
ments for the Goddard Range and Range Rate (GRARR) S-band, modified S-band, very 
high frequency (VHF), and modified VHF systems* The relationships within these four 
systems are analyzed to show that the refraction corrections are the same for all four 
systems and to clarifj^ the group and phase nature of these corrections* The analysis 
is simplified by recognizing that the range rate is equivalent to a carrier phase range 
change measurement. The range errors may be expressed 


AR - ± 


K 



where the + sign pertains to modulation range and the - sign to carrier range, and 
where the equivalent frequencies are 


where 




(carrier) 


f = up link frequency 

f^ = down link frequency 

f - transponder local oscillator frequency 
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ANALYSIS OF IONOSPHERIC REFRACTION ERROR 
CORRECTION FOR GRARR SYSTEMS 

1. INTRODUCTION 

The Goddard Range and Range Rate (GRARR) system is a spacecraft 
tracking system which uses a combination of phase and group ranging techniques. 

Range measurements are accomplished using the basic sidetone ranging 
technique where the phase delay of a ranging sidetone frequency is directly proportional 
to the two-way range between the tracking antenna and the spacecraft. Range rate 
measurements are made using the established principles of coherent Doppler, To 
determine the Doppler frequency, frequency coherence of the ground transmitter car- 
rier frequency is maintained through the spacecraft transponder and back to the ground 
receiver where it is compared in frequency against a coherent sample of the trans- 
mitter frequency. 

The GRARR system exists in two versions, one operating at S-band and one 
operating at VHF. Also, the GRARR stations have been modified so there are essen- 
tially four different systems: S-band, modified S-band, VHF, and modified VHF. 
Ionospheric refraction error correction requirements will be determined for each of 
the four systems. Simplified diagrams of the four systems are shown in Figures 1 
through 4. The diagrams for the S-band and VHF systems were obtained from Ref- 
erences 1 and 3, for the modified S-band and VHF systems, from Reference 2. 

There has been some confusion as to the exact formulation of the ionospheric 
correction for GRARR measurements (Reference 4). This confusion arose because of 
the rather involved nature of the ranging modulation system and uncertainty between 
group and phase treatment. This document is intended to resolve that uncertainty by 
a very careful analysis of basic principles. 

2. PHASE AND GROUP RETARDATION 

In the ionosphere at the frequencies of interest, it is an adequate approx- 
imation to represent the ionospheric range error as inversely proportional to the 
square of the frequency. Thus, a pure sinusoidal or carrier transmitted signal may 
be written in complex form as 

6^)= Re { exp j2rrft] } (1) 
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where 

f = frequency (Hz) 
t = propagation time 

After traveling a distance R, which includes some amount of ionosphere, 
the received signal may be written 

e 2 (t) = Re | exp [j2fff(t - r p )] f 

where approximately for the ionosphere at frequencies above about 100 MHz 

= i J^r _ J|-j = phase transit delay (see Reference 4) (3) 


where 


c = vacuum speed of light 

K is proportional to the total electron content along the path (Reference 5) 

K = 40. 3 f n ds , [meter (Hz)^l (4) 

path 

Q 

n = electron density, e/m 

G 


If the transmitted signal is modulated by any kind of relatively narrowband 
signal, m(t), it can be represented as 

e^t) = Re { exp(j2ffft) exp m(t){ (5) 

where m(t) may, in general, be complex to represent amplitude modulation, phase 
modulation, or both. Then, as shown in Reference 5, irrespective of the form or 
nature of the modulation, the received signal may be written 


where 


e 2 (t) = Re {exp [j2»rf (t - r p )] exp [m(t- r g )] } 
r g - group transit delay 



(6) 

(?) 


For the ionosphere in particular, with r 


as given by equation 3 , 


r 

g 



( 8 ) 


which differs from fp only in the sign of the ionospheric term. The relations of 
equations 3,7, and 8 are used extensively in the following paragraphs to derive 
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the exact ionospheric error dependence of the modulation and carrier phase measure- 
ments for each of the four GRARR systems. 

3. DERIVATION OF SYSTEM EQUATIONS 

For brevity in the followings we adopt a useful shorthand. Phases are ex- 
pressed in millions of cycles (Me) so that the phase of a 10-MHz oscillator iB denoted 
by lot. A phase modulated wave such as 



which is written simply as 

carrier modulation 



modulation-on- 

modulation 

FT 

s 


In accordance with the foregoing conventions , group and phase transit time delays are 



where f , = up and down carrier frequencies 

K - ionospheric range error (positive) normalized to unit frequency. 

Thus, all carrier terms suffer transit delay P u or P d and all modulation or modulation- 
on-modulation terms suffer transit delay or G^ , as appropriate. 
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The derivation for the S-band system is diagramed in Figure 1. The fol- 
lowing notes are keyed to the signal equations at designated points in this figure and 
to the shorthand notation of these equations as presented in Table 1. 

(T) This is the shorthand notation for the transmitted signal, 

A sin(2fff t + fisin2fff t) 

' u r s ' 


© 

© 

© 


The signal arrives at the satellite with time delays and G on the 
carrier and modulation signals. 

The received signal is converted to a low IF by heterodyning with the 

oscillator, Mf . 

o 

The entire resulting low IF signal is used to phase modulate a new 
carrier, Nf Q . Thus the old, carrier ©) becomes modulation and 
the old modulation becomes modulation on modulation. 


© 

© 


The signal arrives at the ground with the down link delays Pj and G_,. 

d d 

This is the N x VCO of a phase-lock-loop, locked to an effective 
60-MHz reference. As such, it does not track any of the modulation 


but insures that the next IF signal will be 60 MHz. By tracing around 
the phase-locked loop, it is seen that <$ 6 must satisfy (denoting 
average or carrier phases by the overbar): 


$ 5 - $ 6 - 50t = lOt 

or <J> 6 = $ 5 - 60t = Nf Q <t - P d ) - 60t 
and subsequently must satisfy 


= $ 5 — $ 5 + 60t = 60t 


© From©, the carrier phase of © is just 60t, the modulation is 
of course, unchanged. 

® *8 ■ V N 


= ('f* 5 - 60t)/N 


(f o N )fc f o P d 


© After this phase demodulator, the old modulation (T) becomes the 
carrier and the old mod -on-mod is the new modulation. 
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(10) This is another phase-loek-loop serving simply to translate the car- 
rier at ( 9 ) to exactly 10 MHz for phase demodulation, ($ 10 = *9 + 
lOt). 

(11) This is the final group demodulation which is phase compared with 
the original signal or modulating frequency f g to give(l2) 

(12) the phase of the group signal (used for the range measurement). 


Meanwhile, 

(^ 3 ) derives the appropriate carrier bias frequency with 
<j> 13 = - (^ - 4> g ) M + 90. 5t, so that at 

( 14 ) we have the carrier Doppler term on a bias frequency, f^, which has 
been derived solely from the ground receiver frequency synthesizer, 
that is, the bias frequency reference is available. (^14 = $13 — ^lO^* 


Equations 12 and 14 constitute the group and phase (or modulation 

and carrier) measurement signals for extraction of range and range rate 


The analysis of the other cases is straightforward and follows a similar 
pattern. 
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4 . 


RANGE AND RANGE RATE EXTRACTION AND EQUIVALENT 
FREQUENCIES 


The modulation phase measurement (My is used for range extraction. In 
each of the four GRARR systems M m in millions of cycles is found to be of the form 


where 


M = 
m 

f = 


G = 

u 


G J = 


f (G +G.) 
s' u d' 

modulating signal frequency 

J7 [ R+ jr] = U P group transit time 
u 

O [ R+ f] - down group transit time 


( 9 ) 


f , f . 
u’ d 


- up and down carrier frequencies 


K = 


thus 


M 


ionospheric range error normalized to unit frequency, (K >0) 
2f_ f K P ^ „ 2f 


m 


s <r r i i i 2 f 

- 5 s ** 

U d 


2f s K 

m 


( 10 ) 


where f ffi is the equivalent frequency for modulation given by 

7 ~ 


( 11 ) 


m 


The form of the apparent range equation is given by solving equation 10 for 
range in the absence of ionosphere errors (K=0): 


m 


2f s, m 


Then the apparent range for K ^ 0 is: 


R 


m 


= R + 


K 


m 


Eq uation 12 is solved t o provide the corrected range measurement 



(12) 


(13) 


Similarly, the carrier phase measurements (M £ ) are used for range rate 
extraction and are all of the form L 

M c - f B* - f L (P d - <V - f u<V G d> (14) 

f B = Bias frequency, see Tables 1 through 4 
= c [ ] ~ up phase transit time 


where 
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p, = -i- R- = down phase transit time 
d c |_ f 2 J 


f = Mf = transponder first local oscillator frequency 
L o 


therefore 


M 


2f Kf F . t 2f 

= f t — R + -^ [4- - 4- + — xl 

B c e Lf 2 f 2 f f 2 -! 

u d u d 

2f Kf r , . f T - f 

= f t — — r + — ~ T4 + 4 + 4 — 51 4] 

bo ' Lf 2 f2 f u n J 

u d a 


-y-StH+Ss. 

13 C .2 


where f is the equivalent frequency for the carrier given by 
c 

_i_ _ i r ± + ± + f L ~ f u ji 

t 2 *Li* i* wJ 


Note thatlf - f |« f and to this extent 
| J-j U| u 

r‘ H? •?] 


(15) 


(16) 


(17) 


As above, the apparent carrier range change equation is the solution of equation 15 
when K = 0: 


* 

R 


c 



u 


M ) 
c' 


when K^O, using equation 15 


* 




(observe the minus sign and compare 
with equation 12) 


( 18 ) 


Equation 18 is solved to provide the corrected carrier range change meas- 


urement. 


R 




(19) 


In principle the corrected range rate measurement is just the rate of 
change of the corrected carrier range change measurement given by equation 19. 
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( 20 ) 


• • * K 

R -R + ~t 
c f 

c 


The actual procedure of range rate extraction is somewhat more compli- 
cated, as described by Zillig, Reference 13. The biased doppler frequency, (M c = 
f b + fj in Zillig's notation, Reference 13) is counted for a preset number, N, of cycles. 
The resulting count period is used to gate the count of a reference or clock frequency, 
f . The resulting clock count, C Q , is the basic measurement recorded. 


C = 
o 


Nf r 


f b + 


■'o ’ 

Nf 


M 


( 21 ) 


where 


* Nf r " C o £ b 

- q is the inferred two-way doppler shift corresponding 


to transmitted frequency, f and includes ionospheric error, 
= bias frequency 
then, from equation 15 


* 

f d 


^ - M c-*b 


2f u A K 

— « E -J2> 

c 


( 22 ) 


This is solved for by the operational equation resulting from setting 
K = 0 , so that for K f 0 we have 


• * c * 

R c = - 2f7 f d = (R 


K 


f 


-) 


(23) 


and the corrected range rate measurement is 

“* K r* * 


„ _ R . K _c * , k _ _ o / N£ r C 0 f b \ 

c t 2 ^' 2f u td 7 c „ — ] 


u 


K 

f 2 

c 


(24) 


In the mechanization of the carrier phase measurement, it is not I ft , as 

c 

implied in equation 22, that is measured, but rather a carrier phase change is meas- 
ured by counting a preset number, N, of Doppler plus bias, f b + f d , frequency cycles. 
The time, t , to count these N cycles is determined by counting the number of cycles, 
C Q , of a higher reference frequency, f occurring during the N-count. 


Thus 




(25) 
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and from equations 15 and 25, 

N * <«b + fj> T * M o <V - M c <*1> 




(R(t 2 )-R(t^ - jr 


* 

As before , this is solved for the measured range change, R, (t 2 ) - 
by setting K = 0 in equation 26, 


R c <*2> - R c <‘l> = 




2f 


V - ( M c 


<y 



Again the asterisked notation is used to denote inferred quantities which include 
spheric error. Then, for K f 0, using equations 15 and 27, 


R* (t 2 ) ' R c (t l> = R (t 2> " R (t l> " 


K (t 2 ) - K (t x ) 


f 2 


c 


Equation 28 is solved to produce the corrected carrier range change. 


R (t 2 > - R (t 1 ) = R c (t 2 ) - R c (t 1 ) 


K (t 2 ) - K (t x ) 


or in terms of measured quantities. 


R (t 2 ) - R (t 1 ) 



+ 


K (t 2 ) - K (t x ) 


f 2 


(26) 

W- 


(27) 

iono- 


(28) 


(29) 


(30) 
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If an expression for range rate rather than range change is desired, it is 
given approximately by dividing the finite range change in equations 29 or 30 by the 
counting interval r . 


R 


R (t 2 ) - R (t x ) 

T 


R*(t 2 ) -R*^) K (t 2 ) - K (t x ) 

+ n 



K(tg)- K (t^) 


f 2 r 

c 


(31) 


that 


This is the operational mechanization of equations 20 or 24, recognizing 


- * 


R 

c 


R c < ( 2> - <y 


T 


and 


♦ 

K 


K (t 2 ) - K (tj) 
— — — 


If there is significant acceleration during the counting interval, so that R 
changes during t and the measured average Doppler frequency f * = - N is in 
general different than the instantaneous Doppler within that interval, then a correction 
term is needed. 


5. RELATION TO PRIOR FORMULATION FOR EQUIVALENT FREQUENCY 

The equivalent frequency formulation that had been in use for operational 
processing of GRARR (Reference 4) carrier data is equivalent to defining f c (old) by 
equation 12 with the positive sign, i.e. 


R* (old) = R + 
c 


K 

(old) 


(32) 
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where f was defined by 
C 1 = 1 
i (Old) ' 1 

The accuracy of the old formulation was questioned when, due to a change in fj and f y , 

equation 33 started to produce imaginary values for f, (old). However, by comparison 

with equation 18, it is not the minus sign of 1/f^ in equation 33 that is in error but the 

plus sign of 1/f 2 . In fact, f (old) should be imaginary if it is defined and used as in 
d c 

equation 32. 

6 , OTHER TRANSPONDER- RECEIVER SYSTEM CONFIGURATIONS 

It is important to note that the results given by equations 15 and 18 for the 
carrier correction are very much a function of the configuration of the transponder- 
receiver system. It happens to be the same (in terms of f u and f d ) for all four of the 
configurations considered here, all of which are basically similar. However, it is 
understood that alternative configurations are under consideration, and a different 
correction may be required for them. 



7 . SIMULTANEOUS USE OF CARRIER AND MODULATION MEASUREMENTS 

TO INFER IONOSPHERIC CORRECTION . 

Due to the dispersive nature of the propagation in the ionosphere (but not 
in the troposphere) it is possible to utilize measurements on two or more frequencies 
in such a way as to eliminate the systematic ionospheric error as in the TRANET Doppler 
system. As an extension of the same principle, it is natural in a system such as GRARR 
to consider the simultaneous use of carrier and modulation measurements to achieve 
the same result. 

From equations 12 and 18, if we compare R m and R c 

R* - R* = KG ( 34 ) 

m c 

where 

G - (l/f| ♦ 1 /£> 
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Since G is known, we can use this relation to infer K, the ionospheric constant through 
the operational equation 


K * = (R m " R c )/G 


(35) 


However, two practical problems arise: that of carrier phase ambiguity 
and that of random errors. In actuality, while the modulation measurement is absolute 
the carrier measurement is only relative, i, e. , subject to an unknown additive con- 
stant. In addition both have random errors that must now be considered. In other 
words, actually 


R* = R 
m 


+ 4 + f 

*2 m 
m 


(36) 


where 


R* = RK+C+f 
c f 2 c 

c 


(37) 


C is a so far unknown additive constant 

f m , < are the random errors having standard deviations a , a . 
me me 


Then solving for K* from equation 35 
K * = K + (f m ‘ C “ ' C )/G 


(38) 


subject to both random ( f and t ) and constant (C) errors. 

me 

Carrier Data Correction . For correcting the carrier data we may now 
take into account that K is used only in differentiated or rate form from which the 
constant C, drops out and is so of no interest* Denoting the approximate differenti- 
ating operator used in the data extraction and preprocessing procedures by D ( ) and 

the result of that operator by the dot notation, i*e. , D (X) = X 


then from equation 18 

• $ 

R c (corn) = D (R * +i£ ) - R* 4 i|— 

c c 


(39) 


and using equation 35 


R c (corr) = 


R ^ 

f?G 


R S 


(f^G - 1) 


(40) 


showing the relative weighting factors for the modulation and carrier measurements. 
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or applying equations (37) and (38) to equation (39) 


R^corr) = R + 


m 


+ <f?Q - 1 )i 


f 2 G 

c 


( 41 ) 


But since to a very good approximation f « f and since there is no particular physi- 
cal reason to expect f and e to be correlated, we have 
r m c 

f 2 G ~ 2 
c 


and 



(42) 


(43) 


Correspondingly, the rms error is 


a^(corr) = 



(44) 


The problem inherent in this is that the random error, < 7 ^, is normally 

very much larger than a • , typically by something approaching the ratio of carrier 

c 

frequency to modulation frequency or bandwidth* In effect, the corrected carrier 
range now has almost ( 1 / 2 ) the error of the modulation measurement, typically two 
to three orders of magnitude worse than the pure carrier measurement. For exam- 
ple,' in the case of GRARR, typical observed short term noise errors are a ^ ~3 
meters and a • « .01 meters/sec (Reference 12) « If, for example, the data were 
differentiated by simple differencing over a one second span, then « Z^2~ = 
4.24 meters/sec. 


This random error degradation can be ameliorated to some extent by 
heavy smoothing on K* as given by equation 35 prior to using it for correction. Clearly 
K* is subject to much heavier smoothing than R* , but the actual improvement available 
would be very much dependent on the actual dynamics of R and K in each particular 
case. If horizontal gradients of ionization level can be ignored, then the constraint 
afforded by the knowledge of the variation of K* with elevation angle could be used to 
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make possible much heavier smoothing, by, in effect, smoothing the vertical column 
integrated refractivity rather than K* which is along the varying slant path. 

Modulation Data Correction . Since it contains an unknown, arbitrarily 
large constant, the estimated ionospheric parameter K* (equation 38) is of no use for 
modulation data correction unless that constant can somehow be estimated. This may 
be feasible assuming the variation of K* with elevation angle to be of known form, then 
curve fitting over the data span yields an estimate of the constant. The full error 
analysis of this regression technique is beyond the scope of this study but will be 
treated in a companion report in this series covering the recovery of ionospheric error 
estimates from dual frequency TRANET Doppler data. 
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Table 1 

S-Band System Equations 

(Phase in Me and Frequency in MHz, See Figure 1) 
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Figure 1. S-Band System, Simplified Diagram 
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Table 2 

VHF System Equations 

(Phase in Me and Frequency in MHz, See Figure 2) 
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Figure 2. VHF System, Simplified Diagram 
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Table 3 

Modified S-Band System Equations 
(Phase in Me and Frequency in MHz , See Figure 3) 



L9 
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Table 4 

Modified VHF System Equations 
(Phase in Me and Frequency in MHz, See Figure 4) 


Modulation 



Modulation On 
Modulation 


f <t - P ) 
u u 

(f - 13f )t - f P 
' u o u u 


12 f o(t- P d) 

(12f + 315)t - 12f P, 
' o o a 


(f - 13f ) (t - G ,) - IP 


f s (t- G u- G d ) 


f (t - G ) 
s' u 

f (t - G ) 

s u' ^ 

(f - 13f )t - f P 
v u O f u u 


<f u -i3f o )(t-a d i-f u p u 
ff u -13f„)(t- G d ) -f u p u f s (t-G, 

A - 13£ c» ft - G d> - ^u,. f B ft - G 1 

£ B <t - G u - G d ) 

£(t-G -G.) 


f (G + G ,) 
s' u d' 


- 13f o + 10)t + 13f o G d - f u< p u + °d 

2 (-6f o + 147.5)t +6f o P d 

3 13f t - 13f P . 

o o a 

4 f u l 

.5 ,(f u - uyt + 13f 0 P d 
6 10t+ 13f (G.-PJ -f(P„ +G d > 


< 10 ‘ 5000>‘ 


it ♦ 131 < G „ - P H> - W + G d> 


= £ B t + f L< G d- p d>- f u ,P u + G d> 
f 

w here f B ■ * ** — 30 kHz , f T ~ i3f 

wnerei B 50()0 * L o 


*Modulation measurement 

** Carrier measurement 
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^TRANSPONDER 11.4166 MhT 
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